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We report on our study of two-flavor full QCD on anisotropic lattices using 0(a)-improved Wilson quarks 
coupled with an RG-improved glue. The bare gauge and quark anisotropics corresponding to the renormalized 
anisotropy ^ = as/at = 2 are determined as functions of /3 and k, using the Wilson loop and the meson dispersion 
relation at several lattice cutoffs and quark masses. 



1. Introduction 

Study of quark gluon plasma using finite- 
temperature lattice QCD has seen a remarkable 
progress over the last decade. The equation of 
state (EOS) is one of the most studied topics, for 
which many calculations have been reported and 
have made clear the effects of dynamical quarks 
jl]. However, most studies on EOS in full QCD 
are made for the temporal lattice size of iVj = 4 
or 6, and a reliable continuum extrapolation of 
EOS is prevented by large scale violation present 
for small Nf. 

Recently, the CP-PACS Collaboration pro- 
posed the use of anisotropic lattice for calcula- 
tions of EOS. This proposal was tested for the 
pure gluon system ||] , for which a well-controlled 
continuum extrapolation of EOS was performed 
for the first time. 

We wish to extend the use of anisotropic lattice 
to calculations of EOS for full QCD. This requires 
a tuning of the bare gauge and quark anisotropy 
parameters which corresponds to a given value of 
the anisotropy. In this paper we report results of 
such a tuning. 



2. Lattice action 

We adopt an 0(a)-improved Wilson quark ac- 
tion coupled with an RG-improved action for 
gluons. This combination of actions has been 
adopted in a series of systematic investigations 
at T = |4j. It also shows better scaling proper- 
ties in finite-temperature QCD, e.g., the expected 
0(4) scaling is observed around Tc 

Here, we extend the study to anisotropic lat- 
tices. We focus on the anisotropy ^ = as/at — 2, 
where as and at are the spatial and temporal 
lattice spacings, respectively. This choice of ^ 
is based on a study of efficiency in calculations 
of EOS as described in our previous study in 
quenched QCD 

We consider the following gauge action which 
includes plaquette, P^^(a;), and rectangle loop, 

Sg = 

Pile' E i^Pvi^^) + + R,^i^)}) 

+jGj2^clPk4{x) + c\Rk4{x) + 4i?4fe(a;))}, (1) 



*Talk presented by T. Umeda. 



where eg -I- 8c^ = 1 and 4 + 4(4 + c^) = 1. 
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Carrying out the RG improvement program of 
Ref. 1^ on this action, we find that a simple choice 
cf = c* = c| — —0.331, which is the same as for 
the isotropic case, sufficiently improves the theory 
for small anisotropy ^ « 1 — 4 [|| . 

The quark action we study is as follows |Q. 



S f 



x.y 



q{x)K{x,y)q{y), 



(2) 



K{x, y) = - nt{{l - n)Ui{x)5,^^i y 
+ (l + 74)(7l(a;-4)<5, . } 



'^s '^{{r - ^i)Ui{x)5 



x+i,y 



i i>j 

The bare anisotropy 7f of the fermion field is de- 
fined by the ratio of the spatial and the temporal 
hopping parameters, Kt — Jf^s- We set the Wil- 
son parameter to be r = 1/^ [Q. The hopping 
parameters are related to the bare quark mass as 
Kg = l/[2(mg 7f + 3r)]. We also define k by 
1/k = 1/ks — 2{'-fF -I- 3r — 4). This parameter 
plays the same role as k on the isotropic lattice. 
For free fermion, equating the rest and the ki- 
netic masses leads to 7^^ = ^^^(1 -I- m^/3) when 
r = l/e@. 

For the field strength i^^;^, we use the stan- 
dard clover-leaf definition. At the tree level, the 
coefficients of the clover terms, Cs and Ct, are 
unity. We incorporate the mean-field improve- 
ment, Ui Ui/uos (i =1,2,3) and — > Un/uQt, 
with the spatial mean-field factor defined with the 
plaquette in 1-loop perturbation theory mqs = 
(1 — 1.154//?)-'^/'*, and Mot — 1 for the temporal 
mean-field. 

3. Calibration procedure 

To realize a consistent anisotropy, the bare 
anisotropy parameters 7^ and 7g for fermion and 
gauge fields have to be calibrated, such that the 
anisotropy S^p calculated in the fermionic sector 
coincides with calculated in the gluonic sector. 

We perform this calibration as follows. Calcu- 
lating (^i?, ^g) for some parameter sets of (7_f, 7g) 



at a fixed n and /3, we make a fit of results with 
linear functions of and 7g, and locate the 
point satisfying ip = £,g = where ^ = 2 in 
our case. This point is denoted as 7!^ and 7^, 
which are functions of (3 and k. 

For the determination of ^g, we adopt 
Klassen's method ^ of matching the Wilson loop 
ratios according to 



Wss{x,y) 



Wss{x + l,y) Wst{x + l,ty 



t = ^GV (4) 



where Wss{x,y) and Wst{x,t) are spatial-spatial 
and spatial-temporal Wilson loop, respectively. 

For ^p we use the relativistic dispersion rela- 
tion of mesons, and demand 



-Q 



(5) 



where E and m are in units of 04, while the spa- 
tial momentum p is in units of Og. The latter 
is defined with p = 2Trn/L,, n = (0, 0, 0), (1, 0, 0) 
and its permutations. We evaluate ^p using either 
pseudoscalar or vector mesons. The results are re- 
ferred to as -f*p{PS) and -f*p{V) (and also 7g(PS') 
and jQiV)) for the calibrated bare anisotropies. 



4. Simulation details 

Simulations are made with the HMC algorithm 
for two fiavors of degenerate quarks, using an 
even-odd preconditioned BiCGStab quark solver. 
The calibration is performed at f3 = 1.8, 1.9 and 
2.0 on a 8^x24, 8^x24 and 10^x30 lattice, respec- 
tively. Six values of k are used at each f3 corre- 
sponding to mps I ray ~ 0.6, 0.7, 0.8, 0.85, 0.9 and 
0.92. Measurements are performed at every 5 tra- 
jectories up to 1000 - 1700 trajectories after 300 
thermalization trajectiries. Errors are estimated 
by the jackknife method with bins of 50 trajecto- 
ries, and error propagation is used in the fitting 
to ^ = 2. 

If the lattice scale is estimated from the Som- 
mer scale rg = 0.5fm, the typical lattice size in 
the spatial direction is about 2.0 fm (2.3 fm) at 

= 1.9, 2.0 (/3 = 1.8) for mps/my ~ 0.8. 
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Figure 1. 1/k dependence of calibration results 
at each (3 

5. Calibration results 

Figure 1 shows the results for and 1/7^^ as 
a function of 1/k at each /3. Unlike the case 
of quenched calculation where l/^p shows 
no linear terms in m^, linear terms are impor- 
tant in full QCD even with the choice r = 
Therefore the results are fitted with a general 
quadratic function of and 1/k, f{(3,K) = 
A + B/P + C/0^ + D/f3K + E/k + F/k^. 

Figure 2 represents the difference between 
jp{PS) and jp{V). Non- vanishing values of this 
quantity represnets an 0{aa) lattice artifacts for 
our action choice. We confirm that the difference 
tends to zero ( less than 5 % for /3 > 2.0) as f3 
is increased. We plan to use this difference to 
estimate systematic errors in continuum extrap- 
olations. 

An extension of the present calibration to the 




Figure 2. Difference between pseudoscalar and 
vector in determination of S,f 

region of weaker coupling (larger /?) is under way. 

This work is supported in part by the Large- 
Scale Numerical Simulation Project of the Science 
Information Processing Center (S.I.P.C) of Uni- 
versity of Tsukuba and by the Grants-in-Aid for 
Scientific Research by the ministry of Education 
(No.12304011, 12640253, 12740133, 13640259, 
13640260, 13135204, 14046202, 11640294). The 
simulations are carried out on Fujitsu VPP5000 
of S.I.P.C. 

REFERENCES 

1. See e.g. S. Ejiri, Nucl. Phys. B(Proc.Suppl.)94 
(2001) 16. 

2. Y. Namekawa et al. (CP-PACS Collabora- 
tion), Phys. Rev. D64 (2001) 074507. 

3. A. Ah Khan et al. (CP-PACS Collaboration), 
Phys. Rev. D63 (2000) 034502. 

4. A. Ah Khan et al. (CP-PACS Collaboration), 
Phys. Rev. D65 (2002) 054505. 

5. Y. Iwasaki, Nucl. Phys. B258 (1985) 141; 
Univ. of Tsukuba Report No. UTHEP-118 
(1983). 

6. S. Ejiri, K. Kanaya, Y. Namekawa and 
T. Umeda, in preparation. 

7. J. Harada et al., Phys. Rev. D64 (2001) 
074501. 

8. T. R. Klassen, Nucl. Phys. B533 (1998) 557. 

9. H. Matsufuru, T. Onogi and T. Umeda Phys. 
Rev. D64 (2001) 114503. 



